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Abstract 

Prompt gluons are an additional source for charmonium suppression in nu- 
clear collisions, in particular for nucleus-nucleus collisions. These gluons 
are radiated as bremsstrahlung in N-N collisions and interact inelastically 
with the charmonium states while the nuclei still overlap. The spectra and 
mean number (ug) of the prompt gluons are calculated perturbatively and 
the inelastic cross section cj^^ is estimated. The integrated cross sections 
a{AB — > J/tpX) for p-A and A-B collisions and the dependence on trans- 
verse energy for S-U and Pb-Pb can be described quantitatively with some 
adjustment of one parameter {ug) cj*^^^. 

I. INTRODUCTION 

The original proposal of using charmonium production in heavy-ion collisions as a 
signal of QGP formation has triggered a series of experimental measurements at the CERN 
SPS in proton- nucleus and nucleus- nucleus reactions(c.f. the reviews [0,0). One remarkable 
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observation out of these experiments is: when one scans through data from the hghter 
to the heavier interacting nuclei, one sees that the J/ip suppression in p-p up to central 
S-U reactions is consistent with a picture of pure final-state absorption on nucleons, while 
the data for central Pb-Pb collisions deviate from this systematics and are therefore termed 
"anomalous" . 

Among the many theoretical concepts we mention the following three mechanisms for ^! 
suppression which have been proposed (here ^ stands for J/'ip^ Xc , i'')'- 

1. Inelastic collisions of a charmonium premeson cc by the nucleons of projectile and 
target. The inelastic cross section o"^b*(''") is a function of the time r (measured in 
the rest system of the premeson until the asymptotic meson is formed) |l7|-[T0|. In 
most calculations an absorption cross section with adjusted value of a^bs = 6 mb, 
independent of time, is used both for the J/ip and ip' channel. While this approach 
reproduces the suppression in p-A collisions and A-B collisions with small projectiles, 
central Pb-Pb collisions show stronger suppressions. 

2. Charmonium dissociation in the quark-gluon plasma phase of the reaction. The inter- 
action of charmonium with the partons, in particular the gluons leads to suppression. 
The different binding energies of the states J/ip, Xc and ip' play a crucial role and may 



lead to discontinuous behaviour in the dependence of the suppression [p|-[T^ 



3. In the final state of a nucleus-nucleus collision many hadrons emerge from the reaction 
volume together with the charmonium. The interaction with these comoving hadrons is 



an additional source of suppression [|T5]-[T^ . In approaches of this kind, the inelastic 



cross section cr^ of the with the comovers is the crucial adjusted parameter for which 
no reliable calculations are available at present. 

In this situation we propose still another mechanism, suppression by prompt gluons. These 
gluons are produced as bremsstrahlung in the N-N collisions in the early phase of the re- 
action. Those gluons, whose production time is sufficiently short, will interact inelastically 



with the charmonium in addition to and at the same time as the \I/-N interactions. In an- 
other language: The gluons produced in a N-N coUision, N + N — > g + N + N may be seen 
together with N as a wounded nucleon A^* = N + g . In this language a \I' will be suppressed 
by collisions with nucleons N in their ground state and by wounded nucleons, A^*. Whatever 
picture one prefers, the prompt gluons will be an additional source of suppression. Their 
number and their momentum distribution can be calculated perturbatively. Uncertainties 
remain for the inelastic cross section o"^^. 

In this paper we expand on the preliminary calculations which were presented in Ref. 
|[T9| . While the general idea is the same, calculations of the suppression by nucleons and 
prompt gluons only, several improvements have been introduced and more observables have 
been calculated: 

(i) A closed expression is derived within a classical multiple scattering model with straight 
line trajectories. 

(ii) Formation time effects are included for the produced premeson and production time 
effects for the gluon by introducing a time dependence for the inelastic cross sections (t) , 
(j^^(t) and for the mean numbers {ng){t) of prompt gluons. 

(ii) We also calculate charmonium suppression in its dependence on transverse energy; 
(iv) We predict suppressions for reactions with inverse kinematics. 



II. FORMALISM 



We follow the original idea of Ref. and derive a closed expressions for charmo- 
nium suppression caused by nucleons and prompt gluons. The derivation is based on Fig.|l|, 
which displays a collision between a projectile nucleus A and a target nucleus B in the two- 
dimensional t — z (time-longitudinal coordinate) plane and in the NN c.m.s.. Since projectile 
and target nucleons have high energies (for the 200 GeV/A energy at SPS, projectile and 
target are characterized by a Lorentz-factor 7 = 10 ), the nucleons move nearly on the 
light-cone. Assuming zero production time for the premeson |T^, a premeson is produced 
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at point 0(^0, -^o) and moves with a velocity Vcm in the c.m.s.. During its motion along a 
straight line trajectory its internal structure develops in time towards the asymptotically 
observed This evolution, called formation, is characterized by a certain formation time 
and manifests itself in a time dependent absorption cross section {t) . It interacts at 
point P{ti,Zi) with a nucleon from the target B. Before this encounter, this nucleon has 
experienced a collision at (5(^2,^2) with another nucleon from the projectile, leading to the 
radiation of a prompt gluon (dashed line) which also interacts with the ^ at P{ti, zi). The 
time difference At^ = ti — to (measured in the NN c.m.s.) has to be put into relation to 
the formation time of the charmonium. Similarly the time difference Atg = ti — ^2 has to be 
compared to the time it takes to produce a gluon (coherence time) , which manifests itself in 
a time dependent mean number of gluons {ng{Atg)). At the point R{t3, Z3) the "if interacts 
with a projectile nucleon Na and its comoving prompt gluons. 

First we treat the interaction of the premeson with target nucleon. The velocity of 
nucleons is taken as 1 (in unit of c). After created at point 0{toZo) but before colliding 
with the target nucleon Nb at point P{ti, Zi), the premeson ^' has lived for a time interval 
Atx^f — ti — to — {z — Zo)/{l + Vcm)- Therefore, the contribution to ^ suppression from the 
target nucleons can be taken into account by an attenuation factor e""^ with 
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FIGURES 




FIG. 1. The collision between a projectile nucleus A with a target nucleus B in the two dimen- 
sional time-longitudinal coordinate representation in the NN cm. system. Nucleons are denoted 
by solid lines, the premeson produced at O(ito,^o) is denoted by heavy dots, while the bremsstrahl 
gluons are drawn by light dashes. 



where we have suppressed the impact parameter dependence of pb, the density of the target 
nucleus B. The time dependent inelastic cross section a^-^ (Ai$) accounts for the formation 
of the charmonium state (see below) . 

We treat the interaction with prompt gluons in the following way: The target nucleon Nb 
collides with a projectile nucleon at point (5(^2,-^2) and a gluon is radiated. This radiation 
process is not instantaneous, but is governed by a characteristic time, the production or 
coherence time (see below). The gluon radiated at (^(^2,2^2) interacts with the \1/ at the 
point P{ti, zi) after a time interval. 




(1) 



Atg = At^Sf 
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The contribution to \I/ suppression from gluons is therefore given by an attenuation factor 
with 

roo /'+00 

Y= dz^psiz)a2{At^) rf/e(Ag7p^(/)air K(Ag) , (2) 

J ZQ J—OO 

where is the inelastic N-N cross section, {ng{Atg)) the mean number of gluons radiated 
in a N-N interaction (to be discussed below), and the 0-function ensures that only gluons, 
which are created before the interaction with the \E', are taken into account. 

The suppression of the \E' via the interaction with projectile nucleons and their accom- 
panying gluons can be derived by appropriate variable exchanges. Combining these results 
together one can write down the \E' suppression factor in A -B collisions as: 

\ C C /'+00 /'+00 

^AB = ^ j d% j d^s j dzAj J dzBj Pa{za,s) pB{zB,b- s) exp[-Ii- I2], (3) 

where Ii is given by 

r+00 

h= dz-f pa{z, s) (Atvp) 

JZA 

*9 /.+00 

X 1 + ^/ dz'^Q{Atg)pB{z\h-s)a^/ {ng{Atg)) , (4) 

and I2 is derived from Ji by the substitutions: za — *■ ^s, Pa{,z, s) pB{z,b — s), Pb{z', b — 
s) Pa{z', s), Vcm -Vcm (hidden in At^). 

The first and second terms in the exponential in Eq. (^ correspond to contributions 
from the nuclei A and B, respectively. In the absence of prompt gluons, (ug) = 0, and 
neglecting the dependence on At^, one has the usual expression for charmonium suppression 
by nucleons. The Et dependence of the charmonium suppression is calculated from the 
impact parameter dependent suppression function 

^ / d^s dzA 7 I-^ dzB 7 Pa{za, s) pB{zB,b- s) exp[-/i - I2] 
J d^s J^^ dzA 7 I-^ dzB 7 Pa{za, s) pb{zb, b - s) 

The suppression for a given value Et is then calculated from 

SUEt) = Jd%P{b;ET)SUb) , 



where the distribution (Pb P(b; Et) gives the probabihty that a given impact parameter 
contributes to a particular value of Ejp. We use the following simplification: The NA50/51 
group has made simulations and give the mean value (P){Et) for each values of Et where 
data exist. We have inserted these values into Eq. @. 

Proton-nucleus collisions can be treated as a special case of Eq. (H). 

For the time dependent effective cross section (At^) between the charmonium state 
^ and a nucleon the expression of the two-channel model for the evolution of a cc pair is 



used 20,10 



(6) 



where afj^ is the asymptotic (r — oo) "^-N cross section, a^^'^ the initial (r 0) premeson- 
N cross section, and AM = M^i — Mjj^. We use values crf/'^^ = 6.7 mb (which includes 



the contributions of Xc and ip' ) and a. 



ip'N 



12 mb. For r ^ we have the absorption cross 



section of the premeson, af^^ , which we take the same in the J/ ip and ip' channels to be 3 
mb lOl. 

The mean number (ug) of gluons radiated in a -A^ collision in the direction of one of 



the two nucleons has been calculated in |jT9[: 



oo fl 
2 



da 



da{qN gX) 
dadk"^ 



(7) 



where 22 



da{qN ^ gX) _ 3as{k^) C 2mla^e + [1 + (1 - af]{k^ + 



dadk"^ 



a + 



91 - a 



a 



Here as{k'^) is the QCD running coupling constant; C is the factor for the dipole approx- 
imation for the cross section of a qq pair with a nucleon, cr^'^(rT) ^ C r^, where is the 



qq transverse separation [^,^,^. C ~ 3 is the pQCD prediction. The coherence time is 
given by 

2Ega{l - a) 



a'^m'^ + k'^ 



and the lower limit of the a-integration 



(t^min + V^^min " ^^)/ (2^q) if k < ^min , 

I k/{2Eq) otherwise. 
Here a is the fraction of the quark light-cone momentum and k the transverse momentum 



carried by the gluon. Compared to Ref. |19| we have introduced a time dependence into 



the mean number of gluons by inserting the square bracket into the integral Eq. (|^). While 
the limits t —>■ and t ^ oo and the characteristic time for the rise are certainly correct, 
the explicit form of the exponential is arbitrary. In our calculation we use k > 0.5 GeV/c 
and thus restrict ourselves to (semi-)hard gluons. Their energy uj > Ci^min has to be chosen 
in relation to the type of charmonium. For example, the destruction of J/tp and ip' needs 
i^min = 0.7 GeV and 0.1 GeV, respectively. However, as long as the premeson are still in 
formation, the uncertainty relation tells us that the binding energies are not yet the final 
ones. 

The choice of tUmin deserves a special discussion since it is related to the problem of 
the energy dependence of "^-g break-up cross section in Eq. (H). At high energies it 



is dominated by gluonic exchange in t-channel and slightly grows with energy pT|. The 
inelastic "^-g cross section is by a color factor 9/4 different from the cross section, which 
is according to the constituent quark model related to the \I/-N cross section. Therefore, one 
expects 

abs ^ abs ^ " abs • 
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0.5 1 1.5 2 

tg (fm/c) 

FIG. 2. Mean number of prompt gluons as a function of the time tg after the NN colhsion for 
various values of Wmin- 



We are left with one adjustable parameter ci^miii in Eq. (0). In Fig. ^we show the mean 
number of gluons, (n^), as a function of length tg for different choices of parameter u^in- In 
the numerical results presented in the next Section we will not use the form of Fig. |^ which 
is generated by the somewhat arbitrary exponential form in Eq. (|^), but use a simplified 
form of {ng{tg)) in order to reduce the numerical efforts: 



{ng{tg)) 



n^gtg/to if tg < to 
if t„ > to . 



n 



This step function type behaviour is also shown in Fig. |^. In our calculations we have used 
a fixed value = 0.6 fm/c, and have varied within certain limits, 0.5 < < 1 for the 
J/ip (corresponding to u^in between 0.5 and 1 GeV) and 1 < n° < 2 for the ip' because of 
its smaller binding energy. 
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III. RESULTS AND DISCUSSIONS 




FIG. 3. The J /tp total cross section as a function of the product AB of the projectile and target 
atomic mass numbers at 200 GeV/c, for = 0, 0.5, 0.75. The data are from p5t| . 

We have calculated the total cross section of J/ip production in proton- nucleus and 
nucleus- nucleus reactions according to Eq. and the results are shown in Fig. |^ for three 
different values of n^. The absolute values of cross section are normalized to the mean value 
of those in p-p and p-d reactions AB= 1 and 2. The long-dashed curve corresponds to the 
result when there are no gluons (only suppression due to nucleons), and cannot reproduce 
the experimental data for the Pb-Pb reactions, while the solid line with = 0.5 accounts 
for all the data very well. 




40 80 120 
Et (GeV) 
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FIG. 4. Ratio of J /ip to Drell-Yan cross section as a function of centrality in S-U collisions 
at 200 GeV/c, with = 0, 0.5, 0.75. The data points are from p. The calculated curves are 
normalized at = to the ratio observed for p-p collisions. 

Various experiments measure charmonium production as a function of the transverse 
energy Et which corresponds to different centralities. Fig.^ shows the calculated suppression 
in S-U collisions as a function of Et with the normalization at Ex = (p-p data). Note 
that the rather large uncertainty introduced by this normalization. 
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FIG. 5. Ratio of J ftp to Drell-Yan cross section as a function of centrality in Pb-Pb collisions 
at 158 GeV/c for n° = 0, 0.5, 0.75. The data are from |^. The calculated curves are normalized 
at Et = to the ratio observed for p-p collisions. 
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FIG. 6. Ratios of J /ip to Drell-Yan cross section as a function of centrality in Pb-Pb collisions 
at 158 GeV/c for n° = 0.5, 0.75 (the two curves in the middle). The upper curve corresponds to 
the result for = 0.5 but normalized at Ex = to p-p data point with a standard deviation 
higher, and the lower one is the = 0.75 curve normalized at Et = to p-p data point with a 
standard deviation lower. 

A comparison of our calculation with NA50 data for Pb-Pb collisions at 158 GeV/c is 
presented in Fig. |^ for = 0, 0.5, 0.75. As for S-U the calculation is normalized at Et = 
to the p-p point. It is obvious that a theory with no gluons cannot reproduce the data. The 
value ?T,° = 0.5 which is predicted for = 0.75 (which corresponds to the J/ip binding 
energy) fits the absolute J/ ip production cross sections (Fig. H) well, but seems a little high 
in the Ej- distribution for Pb-Pb collisions. However, one also has to take into account that 
there is a nearly 15% uncertainty in the normalization (see Fig. 
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FIG. 7. The integrated il^' production cross section as a function of the product AB of the 
projectile and target atomic mass numbers at 200 GeV/c. 
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FIG. 8. The ^ij)' cross section relative to the Dreh-Yan cross section as a function of the transverse 
energy for S-U colhsions at 200 GeV/c @. 



We proceed to the suppression of ip' in nuclear collisions. Since the ip' is less bound 
than the J/ip, we may take oj^iin > 200 GeV and correspondingly inPg ^ 1-2 is the theoretical 
choice and leads to reasonable description of the absolute production cross section (Fig. 
but the Ef distribution in S-U is not well reproduced (Fig.^). The data seem to fall faster 
than the calculation. 

Inverse kinematics 




1 10 10^ 10^ 10'' 10^ 10^ 
AB 



FIG. 9. The integrated J/ij: total cross sections in nuclear collisions in the experimental con- 
figuration {xp = 0.15) and in the situation of inverse kinematics {xp = —0.15) . 
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Et (GeV) 

FIG. 10. Transverse energy distribution of J/ijj production relative to the Drell-Yan for S-U 
collisions {xp = 0.15) and hypothetical U-S collisions {xp = —0.15). 

Finally, we have calculated the integrated and Et dependent cross sections of Jjip pro- 
duction for the inverse kinematics, i.e. p-A and A-p collisions and S-U and U-S collisions. 
We expect differences, since the charmonium is observed with a finite rapidity with respect 
to the NN cm. system. In our calculation we just exchange Vcm —>■ —Vcm of the charmonium. 



The differences arise because of the formation and production times. Figs. ^ and |T0| show 
two results. The effects of inverse kinematics are significant but not dramatic. 



IV. SUMMARY AND CONCLUSIONS 

We have calculated charmonium suppression in nuclear collisions for J/ ip and ip' and for 
integrated cross sections and differential ones with respect to transverse energy. 

The emphasis of the present paper has been on the contribution of the prompt gluons 
to the suppression. This contribution turns out to be significant and decisive, and can 
quantitatively account for the suppressions in almost all cases. However, since the values 
of the input variables n° and cr*^ bear considerable uncertainty, we cannot exclude that 
other mechanisms (plasma, hadronic comovers) will also contribute. In these calculations 
several parameters are adjusted so that the relative importance of the various contributions, 
including those of prompt gluons, cannot be quantitatively determined at the moment. 
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